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Applying a Langevin description of the linear sigma model we investigate four 
different scenarios for the evolution of a disoriented chiral condensate: annealing 
or quench with initial conditions governed by effective 'light' or physical mass 
pions. We present pion number distributions estimated from the zero mode (i.e. 
k = 0-field) component. The best DCC signal is expected for the quench scenario 
with initial conditions centered around (a) fs as would be the case of effective 
light 'pions' close to the phase transition. Our investigations support the idea of 
looking for DCC formation in individual events. 



1 Disoriented Chiral Condensate 

The idea of so called disoriented chiral condensate (DCC) first appeared in a 
workirf AnselmEl but it was made widely known due to BjorkenH and Krzy- 
wickip. This idea has been carried over to the field of high energy heavyrion 
collisions as one of the most interesting suggestions of exotic phenomenaQ as 
it might give direct evidence for the chiral phase transition expected to occur 
at high energy densities. Since then many works appeared oil various aspects 
of DCC formation in heavy-ion collisions (for a review seed and references 
listed). Usually these considerations assume an initial state at high temper- 
ature in which the chiral symmetry is restored by vanishing collective fields 
(or equivalently by vanishing quark condensates). Independent thermal fluc- 
tuations in each isospin direction of the 0(4) sigma-model are present. This 
configuration sits on the top of the barrier of the potential energy at zero 
temperature, so a sudden cooling of the system supposedly brings it into an 
unstable state. This picture is referred to as the quenched approximation □. 
The spontaneous growth and subsequent decay of these configurations would 
give rise to large collective fluctuations in the number of produced neutral pions 
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compared to charged pions, and thus could provide a_mechanism explaining 
a family of peculiar cosmic ray events, the CentaurosQ. A deeper reason for 
these strong fluctuations lies in the fact that all pions are assumed to sit irLthe 
same momentum state and the overall wavefunction can carry no isospinU 

The proposed quench scenario, however, assumes that the effective poten- 
tial governing the evolution of the long wavelength modes immediately turns to 
the classical one at zero temperature. This is a very drastic assumption as the 
soft classical modes completely decouple from the residual thermal fluctuations 
at the chiral phase transition temperature in an ad hoc manner. 

An alternative, the annealing scenario! was suggested by Gavin and Miiller. 
They used the one-loop effective potential instead of the classical one including 
thermal fluctuations. For only moderate expansion in one or more dimension 
this view did not lead to a prediction of a huge DCC signal. The smallness 
of the order parameter at the beginning of the DCC formation has also been 
criticized: if the soft held remains in thermal contact with the fluctuations 
giving rise to the one-loop potential n then one also has to allow for thermal 
fluctuations in the initial conditions BE3. Quenched initial conditions within 
the linear sigma model are statistically unlikely. 

The likeliness of an instability leading potentially to a DCC event during 
the evolution with a continuous contact with the heat-bath of thermal pions 
was investigated by us in til. In the present contribution we will first review 
the theoretical ideas behind a classical Langevin description of the soft fields 
employed inEJ and then shortly discuss our model for simulating the evolution 
of the order parameter and the soft pion fields. We then state some new 
numerical results of the simulation on the evolution of a coherent pionic field 
considered by different scenarios. 

Our picture of a possible DCC formation in high energy heavy-ion collisions 
is as follows: 

• in the first stage of the collision (r = ... 0.5 fm/c) a parton gas is formed 
with a temperature well over the chiral restoration point (T 3> T c ), 

• in the following (r = 1 . . . 2 fm/c) the temperature is around critical (T « T c ) 
and small chirally disoriented domains (bubbles) of collective pionic fields start 
to form together with a thermalized background of (quasi-)pions, 

• at the time (r = 2 ... 10 fm/c) the temperature drops and the fireball expands 
rapidly, at most one bigger DCC domain survives this phase in the heat bath 
of thermal pions, 

• finally (r > 10 fm/c) the fireball breaks off and — in some events — many 
soft pions are emitted with isospin distribution characteristic to DCC. 
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2 Description of soft modes within a thermal environment 

One of the recent topics in quantum field theory at finite temperature or 
near thermal equilibrium concerns the evolution and behavior of the long 
wavelength modes. These modes often lie entirely in the non-perturbative 
regime. Therefore solutions of the classical field equations in Minkowski space 
have been widely used in recent years to describe long-distance properties of 
quantum fields that require a non-perturbative analysis. The distinction be- 
tween soft and hard modes separates the energy and momentum scales as 
gT <C k c <C T, allowing for a perturbative resummation of certain interac- 
tions. Soft modes [k < k c ) are non-perturbative, but - this is the hope - 
can be treated (semi) classically. Hard modes (k > k c ) are treated perturba- 
tively. A justification of the classical treatment of the long-distance dynamics 
of weakly coupled bosonic quantum fields at high temperature is based on the 
observation that the average thermal amplitude of low-momentum modes is 
large and approaches the classical equipartition limit. The classical field equa- 
tions should therefore provide a good approximation for the dynamics of such 
highly occupied modes. However, the thermodynamics of a classical field is 
only defined if an ultraviolet cut-off k c is imposed on the momentum p such 
as a finite lattice spacing a. Many thermodynamical properties of the classical 
field depend strongly on the value of the cut-off parameter k c and diverge in 
the continuum limit (k c — * oo). In a correct semi-classical treatment of the 
soft modes the hard modes thus cannot be neglected, but it should incorporate 
their influence in a consistent way. 

In a recent paper it was shown how to construct an effective semi- 
classical action for describing not only the classical behavior of the long wave- 
length modes below some given cutoff k c , but taking into account also per- 
turbatively the interaction among the soft and hard modes. Also for studying 
non-equilibrium dynamics the separation pf^oft and hard contributions is of 
importance. The resulting effective actionE3E3, S^ff [soft] turns out to be com- 
plex, leading to a stochastic equation of motion for the soft modes. If the hard 
modes are already in thermal equilibrium then the evolution of the soft modes 
is described by a set of generalized Langevin equation - the equation of motion 
corresponding to the above complex effective action. 

We briefly sketch the above following by considering a scalar field with 
interaction £ int = ^j^ 4 - The splitting of the Fourier-components, <j)(p,t) = 
</)(p < k c , t) + <p(p > k c , t), leads to the following interaction part in the action 
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By integrating out the 'influence' of the hard modes on the two-loop level, 
one notices that Feynman graphs contributing at order 0(<? 4 ) to the effective 
action contain imaginary contributions. The stochastic equation of motion has 
the general shape 

□ + m 2 + f^ 3 + J2 f2jV-lV ^" 1(Rer2Jv) ^ = 2 ^ N • (2) 

JV=1 ^ N=l 

Here T2N denotes the effective contribution with 2N soft legs, m 2 the re- 
summed Hartree-Fock self energy (cactus graphs) and £/v are associated noise- 
variables with a correlation (£n £' N ) — ImT 2 jv. When the characteristic time 
scale in the heat-bath of hard modes is short compared to the evolution of the 
soft order parameter fields, then the Markovian limit has a form 

3 

Me//®0 5 + £ (3) 
jV=l,2,3 n=l 

This leads to the familiar form of plasmon damping ('77') and to momentum 
dependent mass corrections. The hard modes act as an environmental heat 
bath. They also guarantee that the soft modes become, on average, thermally 
populated with the same temperature as the heat bath. The noise terms are 
related to the plasmon damping in the Markovian approximation via the simple 
fluctuation dissipation relation 

(£(*)£(*')> = 2TV V S(t - f ) (4) 

in the semiclassical (high temperature) limit. 



3 Langevin description of linear sigma model 

For the description of the evolution of collective pion and sigma fields in the 
0(4)-model we utilize the following simplified equations of motion /taking k c — > 
0) for the order parameters $ a = y J d 3 x 4> a (x,t) in a volume V^lid: 

$i + = &> ( 5 ) 
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with $o = G and $j = (ti, 7T2, 773) being the chiral meson fields and 




(6) 
(7) 



the effective masses. These coupled Langevin equations resemble in its struc- 
ture a Ginzburg-Landau description of phase transition. 

Here "dots" denote the time derivative with respect to proper time r, 77 the 
effective damping coefficient and D/t is the Raleigh-damping corresponding to 
a Z?-dimensional scaling flow of the fireball. The later assumed (rapid) expan- 
sion and cooling results in time dependent temperature and volume according 
to 

T D V D 

1 3t V t 

For the friction coefficient 77 of the a and pion field we take the on-shell 
plasmon damping rate for standard </> 4 -theory generalized to the 0(4)-model: 

where fsp{%) = — fi dtj^j defines the Spence function. We treat the dissipa- 
tion term in a first step as Markovian (mainly because of numerical implemen- 
tation). Furthermore we take m = m n = 140 MeV for specifying the damping 
cofncient r\ in the following. (We remark that both rather crude approxima- 
tions can be relaxed Ej.) In the semiclassical Markovian approximation the 
noise is effectively white and at two-loop level it is Gaussian LL3't3 and thus 
follows a similar relation like (|^). 

Of course, semi-classical Langevin equations may not hold for a strongly 
interacting theory as for highly non trivial dispersion relations the frequen- 
cies of the long wavelength modes are not necessarily much smaller than the 
temperature. Still, when the soft modes become tremendously populated one 
can argue_that the long wavelength modes being coherently amplified behave 
classicallyti Aside from a theoretical justification one can regard the Langevin 
equation as a practical tool to study the effect of thermalization on a sub- 
system, to sample a large set of possible trajectories in the evolution, and to 
address also the question of all thermodynamically possible initial configura- 
tions in a systematic manner. 
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Figure 1: The order parameter <r, the transverse (eff. 'pion') mass my and the longitudinal 
(eff. 'sigma') mass as a function of temperature T for an equilibrium situation for physical 
vacuum pion mass and for vanishing pion mass. 



The 'Brownian' motion of the soft field configuration leads to equiparti- 
tion of the energy at constant temperature. In fig. |l| we show the effective 
transversal masses mj- (^) of the pion modes and tul (R) of the a mode as a 
function of the temperature obtained with solving eqs. (|5j) for a static equilib- 
rium system (i.e. D = 0) at fixed temperature T and sufficiently large volume 
V. The masses shown are thus taken as an ensemble average of the different 
realizations of the Langevin scheme. For large volumes the fluctuations in the 
obtained masses are of the order 1/V and thus small. For the situation that 
the vacuum pion mass is assumed to be zero (no explicite symmetry breaking) 
one can realize from fig. [ljthe situation for a true second order phase transition 
occuring at the transition temperature T = T C = y/2f% — 2ra\j\ w 125 MeV. 
On the other hand for the physical situation of a nonvanishing pion mass of 
TOtt = 140 MeV the 'phase transition' resembles more the form of a smooth 
crossover. In this case, at T w T c , the cr-field still posseses a nonvanishing 
value of (a(T « T c )) » U/2 w a vac /2. 
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4 Stochastic DCC 

In ref. we had studied the average and statistical properties of individual 
solutions of the above Langevin equations with the emphasis on such periods 
of the time evolution when the transverse mass becomes imaginary and 
therefore an exponential growth of unstable fluctations in the collective fields 
might be expected. We had found that in different realistic initial volumes 
individual events lead to sometimes significant growth of fluctuations. 
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Figure 2: The pion number distribution in four different DCC scenarios. 

Here we want to show now the distribution in the number of produced 
long wavelength pions out of the evolving chiral pion order fields within 
the DCC domain V(t) E-i Within the semiclassical interpretation these can 
approximately be written as 

N~° dc « (n 2 + -\ir 2 ) V(r) . (10) 

As a rough estimate one can think of the quench energy density jf^ being 
transformed into low momentum pions yielding N n /V ~ \^ ~ 0.3 /m -3 . 
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On the other hand one might also estimate that the chiral order field 'circles' 
around with (tt 2 + -^r^) < fl resulting in N^/V « / 2 /(2m ;r ) w 0.08 /m" 3 . 
For the total pion number the crucial question is then how large has the evolv- 
ing volume V(t) of the DCC domain increased when the pion oscillations have 
emerged after the roll down period. 

In the following we present results for four different DCC scenarios: 

• Normal or annealing scenario Q : The system is preheated at T = T c for 
sufficiently long time in order to cover a complete set of possible initial thermal 
conditions using its own dynamics (^) and then switch over to a D-dimensional 
scaling expansion. 

• Quench scenario: As for Normal, however, after preheating and thus during 
the expansion the term X/2T 2 of the potential in the equations of motion 
(||,^]) is being omitted in order to mimique an abrupt occurence of the zero 
temperature vacuum potential. 

• m^/f(T = T c ) ks scenario: As for Normal, however, the preheating at 
T = T c is carried out with = (comp. fig. 1) in order to generate initial 
conditions which are centered around the origin a — ir = 0. 

• Quench and m e J^ (T = T c ) w scenario: Like the previous, however, the 
A/2T 2 term of the potential in the equations of motion (||) is being omitted. 

One might argue that the use of the initial conditions prepared within the 
m e J5 (T — T c ) w scenario is inconsistent within the linear sigma model with 
a physical pion mass. From fig. 1 one notices that then the phase transition 
resembles a smooth crossover. From lattice calculations, however, one knows 
that the chiral transitions happens much sharper within a very narrow window 
close at T = T c . This means that it might very well be that the order parameter 
will fluctuate around zero slightly above the critical temperature. 

In fig. 2 we depict the distribution of produced pions within 10f events 
within the four different DCC scenarios. As parameters we had takem D = 3, 
To = 7/m/c and V(tq) = 100 /m 3 . A comparison of annealing and quench 
scenarios, both with finite and vanishing pion mass (for generating the initial 
conditions) reveals that the most productive DCC events would lead to a few 
(6-8 in annealing scenario with finite pion mass), to a moderate number (20 - 
40 in annealing scenario with massless pions and quench with massive pions) 
or to about 120 pions (in quench scenario with initial conditions generated by 
massless pions), respectively. The final distribution does not follow a usual 
Poissonian distribution. Instead, fluctuations with a large number of produced 
pions are still likely with some finite probability! In principle, an ensemble 
averaged description of potential DCC formation carried out within the mean 
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field approximation, as presented in the various literature, can not account 
for such fluctuations and thus has to fail at some point. Also the isospin ratio 
signal in the latter three scenarios is close to that of expected for a DCC event. 
These findings support the idea of looking for DCC formation in individual 
events. 
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Figure 3: The trajectories of the three soft pion modes 7T 1 in the it — cr-plane for each of 
the particular events with the largest pion number out of fig. 2 in the four different DCC 
scenarios. 



In fig. 3 we depict the individual trajectories of the best candidates. Es- 
pecially for the first scenario one recognizes the fact that the initial conditions 
for the order parameters are centered around a « /tt/2 (comp. with fig. 1) and 
do not allow for initial fluctuations in the backward hemisphere with a < 0. In 
any case, all these trajectories clearly illustrate that the collective pion field in 
the late stage of the evolution within the scaling expansion oscillates around 
the 'chiral' circle a 2 + tt 2 ~ f 2 (f n = 93 MeV) as expected for DCC formation. 
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5 Summary and Conclusion 



In summary we want to remark that within the presented Langevin description 
one can simulate, on an event by event analysis, the possible evolution of 
various DCC scenarios in a rather transparent form. It has been demonstrated 
that although the average evolution of the order parameter field is damped, 
there are individual events which show pionic instabilities and appreciable 
pion yields. One can repeat the presented calculations for different sets of 
parameters (the expansion dimension D, the scaling parameter of inital time 
Toli3 and the initial volume V(to). In addition, one can also start to describe the 
evolution with initial conditions generated at much higer initial temperatures 
T T c as carried out by RandrupEj. We will leave this for future workEJ. 
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